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Abstract 

We apply the scale factor duality transformations introduced in the 
context of the effective string theory to the anisotropic Bianchi-type 
models. We find dual models for all the Bianchi-types [except for 
types VIII and IX] and construct for each of them its explicit form 
starting from the exact original solution of the field equations. It is 
emphasized that the dual Bianchi class B models require the loss of 
the initial homogeneity symmetry of the dilatonic scalar field. 
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1 Introduction 



String theory has recently motivated the study of cosmological models 
because its application to cosmology provides an alternative solution to the 
inflationary paradigm, the pre-Big Bang scenario 0. Indeed, the low-energy 
string action possesses a symmetry property, called scale factor duality, which 
lets us expect that the present phase of the Universe evolution is preceded in 
time by a "naturally" inflationary pre-Big Bang phase. In such a scenario, 
the Big Bang represents only the peculiar instant of the Universe evolution 
in which its curvature and its density are maximal. Unfortunately, as far as 
we know, it is not already possible in this context to avoid this maximum to 
be infinite and so, to remove the initial singularity present in the standard 
model. 

Several reasons lead us to think that the standard model characterized by 
an isotropic and spatially homogeneous spacetime cannot be extended until 
the first phases of the Universe 0. Indeed, most cosmologists think that the 
primordial Universe was not necessarily isotropic: for them, it seems more 
natural and more general to say that the Universe began in a less symmetric 
state and became isotropic after some time. So, more general spacetimes than 
FLRW models are often invoked to describe the real dynamical behavior 
of the very early Universe. The cosmological spatially homogeneous and 
anisotropic spacetimes are the most symmetric models after those of FLRW 
type. They are called the Bianchi-type models. 

Since the pre-Big Bang scenario concerns principally the early Universe 
that has to be described by an anisotropic spacetime and since the scale factor 
duality is at the root of this scenario, it is important to demonstrate that it 
remains also valid for any anisotropic spacetime. The aim of this paper is to 
inquire about the possibility of building explicit exact dual solutions - in the 
sense of scale factor duality - for different Bianchi models in order to check 
the validity of the pre-Big Bang scenario in this more general cosmological 
framework. 

The problem of the scale factor duality has already been considered within 
the context of spatially homogeneous vacuum Brans-Dicke cosmologies by 
Clancy et al. 0. This study is based on the Lagrangian formulation of the 
field equations and is restricted to Bianchi type A cosmological models, since 
it is well known that such formulation for class B models is ambiguous M. 
Moreover, the dilaton is supposed constant on the surfaces of homogeneity 
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and the dual and the original equations are assumed to have the same isom- 
etry group, so that, in ||, duality symmetries are found to exist for type J, 
II, VI- X and VH Q . 

In the present paper, we adress this problem in the framework of effective 
string theory^, in a more direct way. We rely on the construction of the dual 
metric from the original one, using the explicit scale factor duality trans- 
formations as deduced from the study of 0(d,d)-invariance of the effective 
action of string theory (JT|, 0)- In particular, we give explicit expressions 
of the dual counterparts of the exact solutions known for different Bianchi 
models in presence of a dilatonic field ||. 

With the exception of Bianchi class A types VIII and IX models, we are 
able to construct explicit dual solutions for each Bianchi spacetime. However, 
for class B models (///, IV, V, VI^-i and Vllh^o), we have to get rid of 
the hypothesis of spatial homogeneity of the dilatonic field and permit its 
effective inhomogeneity character. 

2 The SFD symmetry of the effective action 
of string theory 

In the four dimensional spacetime context, the low energy string effective 
action (in the string frame) can be written asR 



- 2 ( e~+ \R + V a V> - i- H aPl H a ^ 
k z J 12 



^ d 4 x (1) 



Our application of the scale factor duality (SFD) to the Bianchi-type models 
will be made assuming a vanishing antisymmetric field strength H a ^. Within 
this assumption and by varying action ([]]) with respect to the metric g a p and 
to the dilatonic scalar field <p, we find respectively the following field equations 

R a p + g aS V s Vp<j>=0 (2) 

and 

i? + 2D0- [V0] 2 = (3) 

2 Effective string theory with H = corresponds to a special case of Brans-Dicke cos- 
mologies with u = — 1, where to is the Brans-Dicke parameter. 
3 Greek indices always run from to 3. 
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where □ stands for the dalembertian operator. 

It is well known that action (0) is invariant under a SFD transformation 
[[!]]. We have showed recently || that SFD has the same form for any kind 
of metric, i.e. 

G -> G = G- 1 (4) 
0^0 = 0- ln(det G) (5) 

The only difference between a FLRW metric, an anisotropic metric or an 
inhomogeneous metric appears in the building of the matrix G: it always 
contains the metric components relative to the coordinates the metric does 
not depend on, the other components remaining unchanged after duality 
transformations. We can express the only condition necessary for the ap- 
plication of these transformations as follows: "If the metric depends on a 
particular x a coordinate, then we must have 



9af3 







(6) 



for all (3 7^ a" ||. For example, if the metric depends explicitly on x° and 
x 1 and does not depend on x 2 and x 3 , then we can use the transformations 
(Hf) and (H) only if the spacetime metric can take the following form: 
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9oo 
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where G is the following 2x2 matrix: 



G 



922 923 
932 933 



3 SFD in Bianchi-type models 

Bianchi-type models are anisotropic and spatially-homogeneous models 
in which a three-dimensional Lie group of isometries acts simply-transitively 
on the hypersurfaces of homogeneity (for an introduction to the anisotropic 
cosmologies, see e.g. [§). They are nine in number but their classification 
permits to split them into two classes^: there are six models in the class A (J, 

4 In what follows, we shall use the Bianchi models classification presented by Ryan and 
Shepley in @. 
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II, VI-x, VII , VIII and /X) and five in the class B (III, IV, V, VI h and 
V/J/Jf]. As we shall see, for some of these models, the dual dilatonic field loses 
the spatial homogeneity symmetry of the initial field. This is due to the fact 
that, in some casesQ, the determinant of the matrix G becomes dependent on 
a spacelike coordinate and so, in view of (j5|), the initial homogeneous dilaton 
is transformed in an inhomogeneous field after SFD. As will be shown later, 
this loss of symmetry does not affect the spacetime metric which still keeps 
its initial symmetry after duality. 

The nine Bianchi-type metrics can all be written as follows if we adopt 
the hypothesis of diagonality of the spatial metric 

ds * = _ 2 + a 2 2 + b * 2 + c 2 ( ~ 3) 2 (7) 

where a, b et c are the scale factors, functions of the timelike coordinate only 
and where the u a (with a = 0, 1, 2, 3) are the Cartan 1-forms characterizing 
the different Bianchi metrics. But, as we want to apply relations introduced 
previously in ||, we will have to write each metric in its natural basis and 
so to develop Cartan 1-forms in terms of natural 1-forms. 

Most of the string cosmological exact solutions given in the literature are 
not written in terms of the proper time t but rather in terms of the logarithmic 
time t related to proper time by the following differential expression: 

dt(r) = o(r) 6(r) c(r) e^ (r) dr (8) 

So, unless explicit mention, we shall write the metric and perform the du- 
ality transformations in terms of the logarithmic time and use the following 
coordinate system: for any metric, we shall take (x°, x 1 , x 2 , x 3 ) = (r,x,y,z). 
Every initial solution will be noted a, b, c and whereas the barred variables, 
a, b, c and 0, will always refer to a dual solution. 

For some of the class A spacetimes, similar results have been found by 
Clancy et al. || in the framework of scalar-tensor theories. In this paper, 
the authors impose to the dilatonic scalar field to remain homogeneous after 
duality and so they are not able to find a dual solution for each Bianchi-type 
model. Our method to find dual solutions appears however simpler and more 
direct than theirs. 

5 The Bianchi-type model called VIh (resp. Vllh) corresponds to a class B model only 
for ft, ^ — 1 (resp. h ^ 0): for h = — 1 (resp. h = 0), it becomes a class A model. 
6 In fact, in all the class B models. 
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We shall present in the next section a table with explicit exact solutions 
and their corresponding dual solutions for Bianchi classes A and B models. 

3.1 Bianchi / model 

The Bianchi / metric describes the simplest spatially homogeneous aniso- 
tropic model: 

ds 2 = -(abcfe-^dr 2 + a 2 dx 2 + b 2 dy 2 + c 2 dz 2 (9) 

As this metric depends on the timelike coordinate only, the matrix G needed 
for the SFD transformation is the following 3x3 matrix: 

(a 2 \ 
G = b 2 (10) 
V c 2 ) 

with determinant (abc) 2 . From (^), we see that the dual metric can be 
obtained by simply inverting the matrix G: 

ds 2 = -{abc) 2 e~ 2i dr 2 + d 2 dx 2 + b 2 dy 2 + c 2 dz 2 (11) 

and using (|5|), the dual dilatonic scalar field can be written as 

4>= <p-2log(abc) = <p + 2log(abc) (12) 

with, for the dual scale factors, 

a = a~ 1 , b = b~ l and c = c™ 1 (13) 

3.2 Bianchi 77 model 

As mentioned above, the SFD transformations given by (|4]) and ((5|) have 
to be applied to a metric written in its natural basis. The Bianchi // metric 
developed in terms of the natural basis 1-forms can be written as: 

ds 2 = -(abcfe^dr 2 + a 2 (dy - xdz) 2 + b 2 dz 2 + c 2 dx 2 

= -(abc) 2 e~ 2 ^dT 2 + c 2 dx 2 + a 2 dy 2 - 2xa 2 dydz + (a 2 x 2 + b 2 )dz 2 

(14) 
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A spacelike coordinate appears now explicitly in the metric despite its spatial 
homogeneity, so it is necessary to take it into account in the construction of 
the matrix G. As the metric depends on the two coordinates, r and x, G is 
the following 2x2 matrix: 

G — ( ° 2 2 ~? a * 2 ) (15) 
y -xa x a z + b J v ' 

with (ab) 2 as determinant. We note that despite the presence of the coordi- 
nate x in G, its determinant remains time-dependent only. In view of (|5|), 
the dual dilatonic field remains spatially homogeneous as the dual metric. 
In order to retrieve the initial Cartan 1-forms 

uj° = abce~^dr 
oj 1 = dy — xdz 
J 1 = dz 
u 3 = dx 

in the dual metric, we have to add to the inversion of the matrix G the 
following variable change: 

x — > —x 
y -> z 
z -> y 

and so, the dual metric takes the following form: 

ds 2 = -(abcfe- 2 * dr 2 + a 2 (dy - xdz) 2 + b 2 dz 2 + c 2 dx 2 (16) 

with the dual scale factors defined by 

a = b~ 1 , b = a^ 1 and c = c (17) 

The relation (|5]) and the definitions (0) enable one to write the dual dilaton 
in terms of the dual scale factors as follows: 

= (f)-2log(ab) = <p + 2log(ab) (18) 



7 



3.3 Bianchi III model 

The Bianchi III metric is a class B spacetime with the following metric: 

ds 2 = -{abcfe^dr 2 + a 2 dx 2 + b 2 dy 2 + c 2 e 2x dz 2 (19) 

As it depends explicitly on two coordinates, r and x, the matrix G is again 
a 2 x 2 matrix: 

G =( o ) < 20 > 

with determinant (be) 2 e 2x . To retrieve the initial Cartan 1-forms, it is again 
necessary to add the transformation given by (^) and ([5]) the following trans- 
formation on x: x — > —x, with the following result for the dual metric: 

ds 2 = -[abc) 2 e~ 2(S> dr 2 + d 2 dx 2 + b 2 dy 2 + c 2 e 2x dz 2 (21) 

with 

a = a, b — b^ 1 and c = c _1 (22) 
and for the dual dilaton 

<f) = <p-2log{bc) +2x (23) 

This is an example of a dilatonic field becoming inhomogeneous after SFD 
transformations due to the dependence of the determinant of G with respect 
to x. 



Using (|22| ) and (p3|), we can write 

abce~^ = dbce-* +2x (24) 

The expressions of the dual metric and the dual dilaton in terms of the dual 
scale factors are finally given by: 

ds 2 = -(dbc) 2 e- 2 ^- 2x ^ dr 2 + d 2 dx 2 + b 2 dy 2 + c 2 e 2x dz 2 (25) 

= <p + 2log{bc) + 2x (26) 

The inhomogeneity of the metric (|2q) is only apparent, since from (|23T) , 
we can see that the exponential term e~ 2 ^~ 2x ) present in the ^oo metric 
component is only r-dependent, so that the dual metric remains spatially 
homogeneous after SFD transformations. 
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3.4 Bianchi IV model 



The Bianchi IV metric belongs to class B models. As the corresponding 
exact solution we are going to present later is known in terms of the proper 
time and is non-diagonal in Cartan's basis (cf. the table in the next section), 
we shall use the following form for its spacetime metric [10] : 

ds 2 = - a a0 + a 2 a l a l + b 2 j cV(T 2 + d 2^2 a 3 + ^2) + a 3 a 3} ^ 

with a, b, c and d, functions of t and where the 1-forms a % (i = 0, 1, 2, 3) are 
defined by 



a 



a 



dt 
dx 



a 2 = e x dz 



a" 



(dy — xdz) 



(28) 



Developing this metric in terms of the natural basis 1-forms, we obtain: 



ds 2 



-dt 2 + a 2 dx 2 + b 2 e~ 2x {dy 2 + [c 2 - 2 x d 2 + x 2 ] dz 2 + 2 (d 



x 



dydz^ 
(29) 

Again we can see that the metric depends on two coordinates, t and x, 
so that we take for G the following 2x2 matrix: 



G = b 2 e 



-2x 



1 d 2 -X 

d 2 — x c 2 — 2xd 2 + x 2 



Its determinant depends again on the coordinate x: det G = b A e 
Inverting G and making the variable change: 

-x 



-Ax /„2 



(30) 



d 4 ). 



x 

y 



z 

y 



we build the following dual metric: 

as = —o <7 + a a <7 +o <c a a — a [o a + a a ) + a cr f 

with 



a = a, 



b = b- 



z 2 -d A )- 1 '\ 



c and d = d 



and with the same 1-forms a 1 (i = 0, 1, 2, 3) as defined in (p8]). Using 
l2T), we can also write the dual dilatonic field as 



(31) 

(32) 
and 



Ax -log b 4 (c 2 -d 4 ) = (P-Ax + Alog(b) 



(33) 
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3.5 Bianchi V model 

The Bianchi V metric developed in terms of its natural basis 1-forms can 
be written as: 

ds 2 = -{abcfe^dr 2 + a 2 dx 2 + b 2 e 2x dy 2 + c 2 e 2x dz 2 (34) 

The corresponding matrix G is seen to be the following two-dimensional 
squared matrix: 

G = ^(o°) < 35 » 

with determinant e 4x (be) 2 . This determinant being x-dependent, it will lead 
to the inhomogeneity of the dual dilatonic field. 

The building of the dual metric needs both the inversion of the matrix G 
and the change: x — * —x. We can then write the dual metric and the dual 
dilaton respectively as follows: 

ds 2 = -(abcfe- 2 ^- Ax) dr 2 + a 2 dy 2 + b 2 (dz - xdyf + c 2 dx 2 (36) 

= (f) + Ax + 2log(bc) (37) 

with 

a — a, b — b^ 1 and c = c _1 (38) 
Again the inhomogeneity of the metric above is only apparent. 

3.6 Bianchi Vlh model 

Developed in its natural basis, the Bianchi Vlh metric can be written as 

ds 2 = -(abc) 2 e~ 2 ^dr 2 + a 2 dx 2 + b 2 e 2hx dy 2 + c 2 e 2x dz 2 (39) 

For all h ^ —1, this metric is a Bianchi class B spacetime. The particular 
case of h — — 1 transforms this metric in a Bianchi class A model. For the 
purpose of this paper, it is not necessary to consider those cases separately. 
The metric being independent of y and z, G is the following 2x2 matrix: 

n ( b2e2kX ^ , A n\ 
G =[ c 2 e 2x ^ 
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with e 2x( ~ h+1 ) {be) 2 as determinant. 

The dual metric can be obtained by performing both the inversion of G 
and a transformation on x: x — > — x and can thus be written as 

ds 2 = -{dlcfe- 2 * e^ h+ V dr 2 + d 2 dx 2 + b 2 e 2hx dy 2 + c 2 e 2x dz 2 (41) 

with 

a — a, b — b^ 1 and c = c _1 (42) 
Using (^|), the dual dilatonic scalar field takes the following form 

<j> = <p + 2x(h + l) + 2log(bc) (43) 

Note that the dual dilatonic field remains spatially homogeneous only for the 
special case h — — 1, which corresponds to a class A spacetime. 

3.7 Bianchi Vlh model 

We can write the Bianchi Vllh metric in its natural basis as follows: 
ds 2 = -{abefe^dr 2 + a 2 [(X - kY)dy-Ydzf 
+b 2 [Ydy +{X + kY) dz} 2 + c 2 dx 2 
= -{abefe^dr 2 + c 2 dx 2 + [a 2 (X - kY) 2 + b 2 Y 2 } dy 2 (44) 
+2 [b 2 Y(X + kY) -a 2 Y(X- kY)} dydz 
+ [a 2 Y 2 + b 2 (X + kY) 2 } dz 2 

where 

X(x) = e~ kx cos(qx) 
Y(x) = - e~ hx sin(qx) 

q = W^W (45) 

k = h/2 

For all h ^ 0, this metric belongs to Bianchi class B but for the special case 
h — 0, it becomes a Bianchi class A metric. Again the purpose of this paper 
does not require to consider this particular case separately from the others. 

The explicit presence of the coordinate x in the metric implies G to be 
the following 2x2 matrix: 

/ a 2 (X — kY) 2 + b 2 Y 2 b 2 Y{X + kY) - a 2 Y[X - kY) \ 

G ~[b 2 Y(X + kY)-a 2 Y(X-kY) b 2 (X + kY) 2 + a 2 Y 2 ) 

(46) 
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Using the definitions (|45|), we can write the determinant of G as: det{G) = 
a 2 b 2 e~ ikx . 

The dual metric can be obtained making both the inversion of G and a 
transformation on the constant k: k — > —k, i.e. /i — > — /i []. We can thus 
write 

rfs 2 = -(a6c) 2 e- 2 ^ e~ 8kx dT 2 + a 2 [{X - kY) dy - Ydzf 

+b 2 [Ydy + (X + kY) dz} 2 + c 2 dx 2 ' ' ' ' 



a = a , 6 = 6 cm<i c = c (48) 



with 



Using the relation (g) and the definitions (|48f) , we can build the dual dilatonic 
field as follows: 

= 0- 4A;x - 2 log{ab) = (j) - Akx + 2 %(a6) (49) 

Again the inhomogeneity of the g 00 component in the metric ([TJ) is only 
apparent. On the other hand, the dual dilaton remains spatially homoge- 
neous only for k — 0, i.e. for h = (class A model). Indeed, for all h ^ 
(class B model), we lose the initial spatial homogeneity symmetry of the 
dilatonic field. 



3.8 Bianchi VIII and Bianchi IX models 

These are both class A models respectively given (in terms of the natural 
1-forms) by 0: 

• Bianchi VIII metric: 

ds 2 = -{abcf e~ 24> dr 2 + a 2 [2 x dz + (1 - 2xz) dy} 2 

+b 2 [dx + (x 2 - 1) dz + (x + z - z x 2 ) dy} 2 (50) 
+c 2 [dx + (x 2 + 1) dz + (x — z — z x 2 ) dy} 2 

7 This transformation on h is possible because the constant h, in our conventions g], 
has to be comprised between —2 and 2. 
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• Bianchi IX metric: 



ds 2 = —(abc) 2 e~ 2 * dr 2 + a 2 [cos(x)dy + dz] 2 

+b 2 [—sin(z)dx + sin(x) cos{z)dyf (51) 
+c 2 [cos(z)dx + sin(x) sin(z)dy} 2 

In these cases, we cannot apply SFD transformations as given by (Q) and 
(H) in the building of the dual solutions because the conditions (H) are not 
realized. So the method presented in || does not enable one to build dual 
solutions for these two models. 

Nevertheless, it is possible to determine a dual solution for these models 



using Busher's relations given in |TT[ and which can be written, for B = 0, 
as follows: 

9yy = ~ B yy = 

yyy 

g a p = g*p- 9 -^ B a , = o (52) 

9yy 

n — D R - ® ya 

ijya u '-'ya 

9yy 

where y is, for the both models, the only coordinate the metric does not 
depend on and where a and (3 stand for the coordinates the metric depends 
on, i.e. r, x and z. 

It is not necessary to write explicity the two dual solutions to see clearly 
that these SFD transformations introduce a torsion field B which was absent 
initially but, above all, cancel the non-diagonal components of the metric 
so that we lose the initial symmetry of the metric after this SFD. Thus, for 
these cases, the relevance of the SFD transformations is less evident than for 
the others. 



4 Exact dual solutions for Bianchi models 

We shall present in the following table exact Bianchi-type solutions of 
field equations (0) and with their dual expressions. In the first column, 
we shall note the Bianchi-type. In the second and in the third columns, 
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the initial exact solution and its dual will be respectively displayed. All the 
explicit initial solutions presented in the second column come from Batakis 
and Kehagias paper's ]6[] except for the Bianchi IV solution which has been 
obtained by Harvey and Tsoubelis [10|. After any explicit solution, we also 
give the constraint on the constants present therein. Note that the quantities 
N, Pi and qj (with i— 1,2,3 and j=l,2) appearing in the table are constants. 
Introducing both solutions and dual solutions in the field equations (§) and 
(H), we have checked that they satisfy exactly these equations. 

The Bianchi Vllh as well as VIIq general exact solutions of field equations 
(0) and (|3p being not known, we have not been able to build their dual 
counterparts in explicit form. 
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Type 


Solution 


Dual solution 


I 


a(r) 2 = e { P^ T 

6( r )2 = e fe+JV)r 

c (r) 2 = e ( p3+JV ) r 
0(r) = iVr 


a(r) 2 = 
Krf = 
c(r) 2 = 
0(r) = 


e -(pi+N)r 
6 ~(P2+N)t 

-(2N + Pl + p 2 + p 3 ) r 


with 


Ei<j ft Pi = N 2 






II 


a(rf = X- l e NT 
b{rf = X e^ +N)T 
c(r) 2 =Xe^ 2+N ^ 
0(r) = Nt 


a(r) 2 = 
b(r) 2 = 
c(r) 2 = 
0(r) = 


X' 1 g-(2pi+AT)r 
^ e (2p 2 +A r )r 

-(AT + 2pi)r 


with 


X(t) = l/p 3 cosh(p 3 r) 


4»iP2 — 


p 2 = TV 2 


III 


a(r) 2 = Pl e fe+JV)r sinh -2( pi r ) 

6(r) 2 = Pl e (7V ~ P2)r 

c (r) 2 = Pl e fe+^)- sm /i- 2 ( pi r ) 

0(r) = A^r 


a(r) 2 = 
6(r) 2 = 
c(r) 2 = 
0(r) = 


pi e (P2+JV)r sinh -2( pi r ) 
p -l g(P2-JV)r 

p- 1 smh 2 {p l r) e -fe+^)- 
-JVr + 2x -2log(p 1 ) 
+2 /og (sinh(pi r)) 


with 


/I 2 2 7\r2 

4 Pi — P2 = N 






IV 


a(t) 2 = (a'H) 2 
b(t) 2 = (a Hf a0 
c(t) 2 = 1 + a /og (a 1 tj 
d(t) 2 = a log (c^t) 

m = 


a(t) 2 = 
6(t) 2 = 

c{tf = 
d{tf = 

$(t) = 


(«„-'t)" 2w 

1 + a ^ ( a o 1 

— 4x — 4 a log ( a o 1 ^) 


with 


a = 4/5 
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Type 


Solution 


Duale 


V 


y ' 2 

c(r) 2 = 1 Xe^ P2+N ^ T 
y ' 2q 2 

0(r) = iVr 




a( T ) 2 = ^Xe NT 
2 

b(r) 2 = 2 e -^ 2+N > 

c(rf = 2q2 e^~ N > 
X 

4>(t) = -2log(X(r)) - Nt 
+4x + 2log(2) 


with 


X(r) = Pi / sinh(pi r) 




3p 2 -p 2 = N 2 


VI-! 


a(r) 2 = qj Pl exp [v\e 2p2T ] e {pi+N ^ T 
b(r) 2 =p 2 3 p 2 e {p2+N)r 

c(r) 2 =p 2 3 p 2 e (p2+N)T 
0(r) = Nt 


a(r) 2 = q\p x exp \v\t 2p2T \ e^ Pl+N ^ T 
Kr) 2 = {p\p 2 y l e^ P2+N > 

c(r) 2 = {PIP2)' 1 e~^ 2+N > 

0( r ) = -(N + 2p 2 )r-2log(p 2 p 2 ) 


with 


p 2 2 + 2p lP2 = N 2 






VI h 


a(r) 2 = q 2 1 X^e^ P2+N > 
b(r) 2 = q 2 2 X h{h+ ^ e^+^ T 
c(r) 2 = q 2 2 X^ e ( - p2+N ^ T 
0(r) = Nt 


a(r) 2 = g 2 X (ft — P2+A M T 
6(r) 2 = 5 -2 X -Mft+i) e -(M+tf)T 
c( r ) 2 = g 2 X-( ft+1 ) e (^-JV)r 
0(r) = -7V T + 2 x (/i + 1) 

-(/i + l) 2 W[X(r)] 


with 


X(r) = 


smhvpi t) 

Pi 


-2 
(h+1)' 2 


(h+1) 2 1 2 
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